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ABSTRACT 

Recent observations that indicate that some extrasolar planets observed in 
transit can experience mass loss from their surfaces. Motivated by these find- 
ings, this paper considers outflows from Hot Jupiters in the regime where the 
flow is controlled by magnetic fields. Given the mass loss rates estimated from 
current observations — and from theoretical arguments — magnetic fields will 
dominate the flow provided that field strength near the planet is greater than 
~ 1 gauss, comparable to the surface fields of the Sun and Jupiter. The prob- 
lem can be separated into an inner regime, near the planet, where the outflow is 
launched, and an outer regime where the flow follows (primarily) stellar field lines 
and interacts with the stellar wind. This paper concentrates on the flow in the 
inner regime. For a dipole planetary field with a spatially constant background 
contribution, we construct a set of orthogonal coordinates that follow the field 
lines and determine the corresponding differential operators. Under the assump- 
tion of isothermal flow, we analytically find the conditions required for escaping 
material to pass smoothly through the sonic transition, and then estimate the 
mass outflow rates. These magnetically controlled outflows differ significantly 
from previous spherical models: The outflow rates are somewhat smaller, typi- 
cally M ~ 10 9 g/s, and the flow is launched primarily from the polar regions of 
the planet. In addition, if the stellar wind is strong enough, the flow could be 
reversed and the planet could gain mass from the star. 

Subject headings: MHD — planetary systems — planetary systems: formation 
— planets and satellites: formation 
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1. Introduction 

Among the hundreds of extrasolar planets discovered to date, a substantial fraction orbit 
their stars with periods of 10 days or less. These planets are thought to have formed further 
out in their solar systems (e.g., Lissauer & Stevenson 2007), and subsequently migrated 
inward (e.g., Papaloizou & Terquem 2006) where they become stranded at small semimajor 
axes, perhaps due to disk truncation (Shu et al. 1994, Lin et al. 1996) or because the disk 
loses so much gas that it can no longer move planets. After planets reach these inner orbits, 
they are subjected to intense heating from their parental stars. This heating, which is most 
effective for UV photons, can drive photo-evaporative flows from the planetary surfaces. 
In the most extreme cases, the resulting mass loss could affect both the final masses and 
densities of the planets. In other cases, the outflows can observable — even if their effect 
on the final mass is modest - - and can provide important information about planetary 
properties. This paper explores outflows from Hot Jupiters in the regime where magnetic 
fields are strong enough to guide the flow and thereby determine the outflow geometry. 

Mass loss from Hot Jupiters has been observed in association with the transiting planet 
HD209458b (Vidal-Majar et al. 2003, 2004; see also Desert et al. 2008, Sing et al. 2008, 
Lecavelier des Etangs et al. 2008). A recent followup observation, using the Cosmic Origins 
Spectrograph on the Hubble Space Telescope implies a mass outflow rate M^8x 10 10 g/s 
from this planet (Linsky et al. 2010). In addition, signatures of atmospheric evaporation 
from the extrasolar planet HD 189733b have recently been reported (Lecavelier des Etangs 
et al. 2010). More examples are expected in the term future, as well as null detections 
(Lecavelier des Etangs 2007). As a result, the collection of close-in extrasolar planets provides 
a laboratory to study the process of mass loss from planetary bodies. 

Theoretical models of mass loss from extrasolar planets have been considered previously. 
Pioneering models of outflows from these planetary bodies have been constructed (Lammer 
et al. 2003; Baraffe et al. 2004, 2006) and indicate that substantial mass loss can take place. 
However, related studies of the effects of mass loss on the population of close-in extraso- 
lar planets show that it is difficult to explain the observed mass distribution (Hubbard et 
al. 2007). In any case, a number of open questions remain. The aforementioned studies 
(primarily) use simple energy-limited outflow models (see also Watson et al. 1981), in con- 
junction with physically motivated scaling laws, and produce a range of outflow rates for 
given planetary masses and external UV fluxes. The next generation of theoretical calcu- 
lations considered refined treatments of the chemistry, photoionization, and recombination 
(Yelle 2004; Garcia Muhoz 2007), as well as the effects of tidal enhancement (Murray-Clay et 
al. 2009, hereafter MCM). More recently, two-dimensional models of planetary winds have 
been considered (Stone & Proga 2009) and indicate that the mass loss rates can be less than 
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those in the spherical limit. In addition to planetary outflows, alternative explanations of 
the observations have been put forth, where the inferred excess material is due to a confined 
exosphere (Trammell et al. 2010; hereafter TAL) or a mass transfer stream (Lai et al. 2010). 

This theoretical investigation into planetary outflows adopts a new approach. As out- 
lined below, the magnetic fields - from both the star and the planet - are generally strong 
enough to guide the flow (see also TAL). When the outflow follows the magnetic fields lines, 
the geometry of the flow pattern is set by the field structure, but can be quite complicated. In 
particular, the outflows depart significantly from spherical symmetry and existing (primarily 
spherical) models are not applicable. In spite of this complication, the outflow problem can 
be reduced to one dimension by constructing a new orthogonal coordinate system where one 
coordinate follows the magnetic field lines. This approach allows for the outflow properties to 
be determined semi-analytically. In this context, the term "semi-analytic" refers to models 
where the equations are reduced to, at most, ordinary differential equations. 

The problem of magnetically controlled outflows from planets can be divided into sub- 
problems. Since the magnetic field structure of the planet guides the flow — for the regime 
considered here — the field is independent of the outflow and can be determined separately. 
Section [2] discusses the conditions required for this approach to be valid and provides an 
overview of the relevant scales in the problem. This work assumes that the planetary field 
has a dipole form. Near the planetary surface, the stellar component of the magnetic field is 
smaller than the planetary field and is slowly varying (spatially); as a starting approximation, 
we thus assume that the stellar field provides a constant background contrribution. For this 
geometry, we construct a set of orthogonal coordinates that follow the field lines and hence 
the flow (Section [3]). For a given field configuration, the outflow problem can be separated 
further into two regimes: [1] The launch of the outflow near the planetary surface; this flow 
depends on the heating and cooling in the vicinity of the planet, and the flow geometry is 
constrained by the magnetic field configuration. [2] The propagation of the outflow, away 
from the planet, as enters the regime where the magnetic field structure and gravitational 
forces are dominated by the star. This paper focuses on the launch of the wind, in Section HJ 
and the propagation problem is left for future work. The observational implications of these 
outflows are discussed in Section The paper concludes, in Section [HI with a summary of 
results, a discussion of their implications, and some directions for future work. 

2. Overview 

This section defines the basic scales in the problem and justifies our approach. After 
presenting an estimate for the outflow rate, we show that the gas is well-coupled to the 
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magnetic field and that the magnetic pressure is larger than the ram pressure of the flow 
by several orders of magnitude. We also discuss the background magnetic field provided by 
the star and the effects of the stellar wind, where both of these effects limit the sphere of 
influence of the planet. A background magnetic field component (from the star) is included, 
whereas additional contributions from currents are shown to be small. 

We first make an order of magnitude estimate for the mass outflow rates from Hot 
J up iters using a simple scaling argument: If we assume that the outflow is limited by the 
rate at which the gas gains energy from the stellar UV flux, then the mechanical luminosity 
of the outflow GMpM / Rp must balance the rate of energy deposition, 7] ia dFuv^Rp, where 
the parameter ?7 ra d includes the efficiency of energy capture and takes into account the fact 
that radiation can be absorbed above the planetary surface (at Rp). The resulting mass 
outflow rate M is thus given by 

where this expression would be exact if one could determine the correct value of ?7 ra d. Next 
we note the coincidence that a particle of mass fi living at the surface of a Jovian planet has 
potential energy e M given by 

f „ = ^«13.9eV, (2) 

Up 

where we have used typical values /i = m P) M P = lMj, and R P = 10 10 cm to evaluate the 
energy. This potential energy scale is almost the same as the ionization energy for Hydrogen 
atoms. As a result, we can write the mass outflow rate in the form 



450 erg s -1 cm -2 



(3) 



where the first (dimensionless) term in brackets is close to unity and where the second term 
in brackets represents the number of UV photons intercepted by the outflow per unit time. 
The UV flux is scaled to the benchmark value F uv = 450 erg s -1 cm~ 2 , the flux appropriate 
for the quiet Sun at a distance of a = 0.05 AU (e.g., MCM, Woods et al. 1998). The 
resulting numerical estimate for the outflow rate agrees with previous results (e.g., Garcia 
Munoz 2007). When the efficiency r] rad is high, the outflow rate is thus determined by an 
approximate balance with one outgoing particle per incoming UV photon. 

The above estimate ignores magnetic fields, whereas this paper considers the wind to 
be guided by the field. In order for the plasma to be well-coupled to the magnetic field, 
the cyclotron frequency uc must be larger than the collision frequency T. The cyclotron 
frequency is given by uc = Bq / (mc) and the collision rate is given by T = nuv. For 
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the parameter space of interest, we expect the magnetic field strength near the planet to 
be B ~ 1 — 40 G (e.g., TAL) and the collision cross section to be a ~ 2 x 10~ 13 cm 2 
(see Shu 1992, Spitzer 1978, Surrock 1994). Because the continuity equation implies that 
M = 4nr 2 pv, we can write the ratio of frequencies in the form 



u c qB AnqBr 2 A f B\ ( a \~ X I M 



T cmnav ca M " W \IG ) \2 x lO" 13 cm 2 ) yiO 10 g s" 1 ' ( l) 

Although the field strength decreases with distance from the planet, equation (jlj) shows 
that the frequency ratio scales like Br 2 . Since the dipole field strength B oc 1/r 3 , the ratio 
wc/r decreases as one power of the radius and hence formally exceeds unity out to a radius 
rc ~ 10 4 i?p, well beyond the launching radius r s of the wind (where r s « 3Rp; see Section 
H]). This estimate for rc is much larger than the radial scale where the stellar environment 
- including the stellar wind, magnetic field, and gravity — dominates that of the planet. 
Although the number density of the wind may continue to decrease in this regime, the 
magnetic field strength will be larger than the scaling used here and the wind will remain 
tied to the field. We thus conclude that u>c 3> T for the regime of parameter of interest, and 
that the outflow is well coupled to the magnetic field. 

Another necessary condition for the magnetic field to guide the outflow is that the 
magnetic pressure must be larger than the ram pressure of the flow. Here we find the radius 
where the two pressures are equal, i.e., where pv 2 ~ B 2 /87r. We can write the density in 
terms of the outflow rate M ~ 4iir 2 pv = T^itr 2 pa s , where a s is the sound speed and where 
the second equality defines the parameter J r m . The magnetic pressure and the ram pressure 
are equal when B 2 r 2 = 2Ma s J r ra . After scaling the magnetic field using the usual dipole 
relationship, so that B = Bp(Rp/r) 3 , we can for solve the radius within which the magnetic 
field is dominant, i.e., 

r ,„ _ \ _1 / 4 „ / Bp R P \ ' M 



(BpRpf 2 (2Ma s F m ) w 26 l-^ 



R P y r r ' \ "7 ViOG 10 10 cmy \l0 10 g/s lOkm/s 

(5) 

This equation is evaluated using the surface field Bp = 10 G, the value estimated for ex- 
trasolar planets (Christensen et al. 2009; see also TAL); for comparison, the surface field 
strength for Jupiter is somewhat lower, Bp = 4.2 G (Stevenson 2003). This result implies 
that the magnetic pressure exceeds the ram pressure of the outflow for radii near the planet, 
where "near" is defined to be within about 26 Rp. However, the stellar magnetic field is 
generally stronger than the planetary field at these radii, so that the magnetic field pressure 
dominates at all distances from the planet. In addition, the sonic transitions take place at 
r ~ 3Rp, well within the boundary defined by equation (jSJ). Using the dipole scaling relation 
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for the field stength, we find that the magnetic pressure is larger than the ram pressure of 
the outflow by a factor of ~ 10 4 at the sonic surface and by a factor of ~ 10 6 at the planetary 
surface. As a result, the outflow must be magnetically controlled. For completeness, we note 
that these field configurations contain an X-point, a location where the field vanishes in the 
equatorial plane and where the condition ((Sj) for magnetic pressure domination fails; this 
complication affects only a few streamlines and reduces the overall mass outflow rate by a 
small amount. 

The star supports a stellar wind that also provides a ram pressure. In order for the 
planet to successfully launch an outflow under magnetically controlled conditions, the ram 
pressure of the stellar wind cannot be too large. The ratio of the ram pressure of the stellar 
wind to that of the planetary outflow can be written in the form 



where r is the radial coordinate centered on the planet and the stellar parameters are eval- 
uated at the location of the planet (a distance a from the star). If the stellar wind has an 



for the Sun, then the ratio of ram pressure in equation ([6]) is close to unity within a few 
planetary radii. As shown above, B 2 ^> pv 2 for the planetary wind, so that the same is true 
for the stellar wind; as a result, the magnetic pressure near the planet will be much larger 
than the ram pressure of both the stellar wind and the planetary wind, and the outflow will 
be magnetically controlled. At large distances from the planet, as determined from equation 
fl5J, the stellar wind pressure will play an important role. 

The star also has a magnetic field that must be taken into account. If the surface strength 
of the stellar field is B*, the strength at the location of the planet will be B ps B*(R*/a) 3 , 
where a is the semimajor axis of the planetary orbit. This expression would be exact in 
the limiting case of a pure dipole field with the planetary orbit in the equatorial plane 
of the dipole. In practice, however, the stellar field will be more complicated due to the 
stellar wind. For solar- type conditions, the field lines are closed out to radii ~ 3i2*, but 
they spiral outwards at larger distances. The planet radius Rp ~ 10 10 cm is much smaller 
than both the stellar radius R* ~ 10 11 cm and the distance to the star a ~ 10 12 cm. As 
a result, to leading order, the stellar contribution to the field can be considered to have 
constant field strength and constant direction over the region where the planetary outflow 
is launched; it is straightforward to show that this approximation results in error terms of 
order 0{r/a) ~ 10~ 2 . At the planetary surface, the ratio of the stellar field strength to that 
of the planet is given by 




(6) 



outflow rate M* = 10 12 g/s and outflow speed v* 



400 km/s, comparable to the values 




(7) 
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where is the field strength on the stellar surface and B(a) is the stellar field strength 
evaluated at the position of the planet. For the regime of parameter space of interest, this 
field geometry has an X-point, a location where the magnetic field vanishes (Shu 1992). Here, 
the stellar field is (nearly) constant in the vicinity of the planet, whereas the planetary field 
decreases as B oc 1/r 3 . To leading order, the radius of the X-point is given by 

r X ^^(^-X /3 -10R P . (8) 



For comparison, note that the Hill radius is given by r# = a(Mp /SM*) 1 ^ . For a Jovian 
planet and solar-type star, th ~ 0.07a ~ 7Rp. As a result, the sphere of gravitational 
influence of the planet (determined by m) and the sphere of magnetic influence of the 
planet (determined by r x ) are approximately the same. 

To illustrate the magnetic field geometries that arise, we plot a representative collection 
of magnetic field lines in Figure [TJ This plot uses (3 = 0.0014, a value that occurs (for 
example) when the star and planet have comparable surface field strengths and the length 
scales R*/Rp = 10 and a/Rp = 100. Magnetic field lines that originate from small polar 
angles on the planet surface curve off to large distances. These field lines will, in general, 
connect up with stellar field lines, which either end on the star or are carried out to large 
distances by the stellar wind. Field lines that start at lower latitudes (closer to the planetary 
equator) curve back and end on the planetary surface, i.e., they are closed. For a given field 
configuration, a well-defined fraction Fap of the planetary surface is exposed by having its 
field lines open (relative to the planet). The field line that passes through the X-point, shown 
by the dashed curve in Figure [U delineates the boundary between the closed and open field 
lines. 

In principle, currents produced by the flow could modify the planetary magnetic field, 
assumed here to have a nearly dipole form. This effect has been extensively studied, primarily 
in the context of stellar winds (e.g., Mestel 1968, Mestel & Spruit 1987; see also TAL 
for an application to extrasolar planet magnetospheres) . As shown in the aforementioned 
work, currents J± that run perpendicular to the flow can arise from either vorticity or 
from the variations of the Bernoulli constant across field lines. These currents give rise to 
perturbations in the magnetic field strength, B±, which are of order 

B ± Airpv 2 



B ' (9) 

As shown above, this ratio is quite small within the sonic surface, typically less than 10 -4 , 
so that the additional fields produced by perpendicular currents can be neglected. 



Magnetic fields can also be generated by current sheets, which arise at the interface 
between the field lines that carry outflowing material and those that do not. The difference in 
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Fig. 1. — Magnetic field lines centered on the planet for field configurations with background 
parameter /3 ~ 0.0014. The solid curves show a set of field lines for starting values of sin# 
that are evenly spaced. The dashed curve shows the field line that passes through the X- 
point. The field lines near the poles are open, and join onto the background field lines 
provided by the star and stellar wind. The field lines near the equator are closed and end 
back on the planetary surface. In this coordinate system, distances are measured in units of 
the planet radius, and the star is located at position (x,y,z) = (100,0,0). 
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flow velocity across this boundary implies a pressure difference (from the Bernoulli equation) 
which results in a perturbation SB of the magnetic field. Applied to the present problem, the 
results of TAL imply that (SB)/ B ~ 87rP/B 2 , where P is the thermal pressure and B is the 
dipole field of the planet. Over most of the parameter space of interest, the field produced 
by current sheets (SB) is thus the same order as that produced by perpendicular currents 
(B±, see equation [9]), and hence can be neglected to leading order ([SB]/ B ~ 10~ 4 ). 

Perpendicular currents and current sheets can be neglected when the pressure P <C 
B 2 /(8ir). In the case of stellar winds, this condition is violated sufficiently far from the 
stellar surface; at these large distances currents cannot be ignored, as they lead to open 
field lines, which are necessary for winds to escape (Mestel & Spruit 1987). In the present 
case, however, an external field allows for open field lines (relative to the planet) even in the 
absence of currents. As shown in Section 3.2, as long as the external field is nonzero, some 
of the field lines will be open; the fraction of the planetary surface that is covered by open 
field lines is an increasing function of the external field strength (see equation [27]). 

The above discussion defines the expected planet radius Rp, X-point radius r x , Hill's 
radius r#, stellar radius i?*, and semimajor axis a. The principal calculation of this paper 
determines the radius r s of the sonic point, the location where the outflow is launched and 
can escape the planet, and finds r s ~ 3 — ARp (Section H]). The remaining length scale 
in the problem is the scale height H of the planetary atmosphere, where H = kT / (fig) = 
kTRpf(fiGMp). As a reference point, we evaluate the scale height at the atmospheric level 
where most of the stellar light is absorbed. At this layer, previous work (e.g., MCM) shows 
that the effective temperature of the planet T ~ 1000 K and hence H ~ Rp/100; the scale 
height increases to H ~ Rp/10 in the upper atmosphere where the UV photons are absorbed 
and T ~ 10 4 K. With this specification, the length scales involved in planetary outflows obey 
the ordering 

H < R P < r s < r H ~ r x ~ R* « a . (10) 

The sphere of gravitational influence of the planet (from the Hill radius r#) and the sphere 
of magnetic influence (from the X-point radius rx) are roughly the same size, and are 
comparable to the stellar radius. The sonic surface typically lies at a few planetary radii and 
is thus well inside both r# and rx- On the other hand, the sphere of influence of the planet 
is much smaller than the star-planet distance (the semimajor axis a). This separation of 
scales allows for the launch of the outflow to be considered independently of the subsequent 
propagation of the flow. 
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3. Magnetic Field Geometry with Constant Background 

To gain further understanding of this problem, and to simplify the calculations, we use 
the method of matched asymptotic expansions. Specifically, we divide the the problem into 
two regimes: [1] The region near the planet where the outflow is launched, and [2] the region 
"far" from the planet where the magnetic field structure is determined by the field of the star 
(modified by the stellar wind). In the near region, the magnetic field is primarily determined 
by the dipole field of the planet, but nonetheless contains a contribution from the stellar 
field. However, this stellar contribution is nearly uniform. Exploiting this property, we can 
model the magnetic field near the planet through the reduced form 

B = Bp [C 3 (3 cos 6r-z) + f3z] , (11) 

where £ = r/Rp and where (3 is defined by equation (j7j). For simplicity we have taken the 
background field to point in the z direction. Notice that this reduced field is axisymmetric, 
so that the problem becomes two dimensional. In addition, the field is anti-symmetric with 
respect to reflections across the z = plane; since the sign of the field does not affect the 
dynamics, the outflow is the same for both hemispheres of the planet. 

In this initial treatment, the dipole field of the planet is augmented by a constant 
background field (due to the star). The stellar field is expected to be dipolar near the 
stellar surface, but will be modified by the stellar wind. This stellar wind acts to straighten 
the field lines, making them more radial. In any case, the stellar field is expected to be 
nearly constant, in both strength and direction, over the region where the wind is launched. 
However, the background field will not necessarily point in the z direction of the planetary 
dipole (as assumed here). Notice also that the planetary spin, and hence the planetary 
dipole, does not necessairly line up with direction of the orbital angular momentum, i.e., the 
dipole also has (in general) an arbitrary direction (e.g., Fabrycky & Winn 2009). This paper 
considers the simplest case where the dipole and the background field are aligned (equation 
fTTJ). Despite its simplicity, this ansatz displays the key features of the expected magnetic 
configuration: a dipole form near the planet, an (effectively) straight geometry far from the 
planet, and some open field lines even in the absence of currents. The more general case, 
where the two components lie at an arbitrary angle, should be considered in future work. In 
addition, if the star rotates more slowly than the angular velocity of the planetary orbit, the 
field lines will tend to wrap up into spiral configurations (Parker 1958); this effect should 
also be included in future studies. 
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3.1. Construction of the Coordinate System 

With this configuration, the magnetic field is current-free and curl-free, and hence can 
be written as the gradient of a scalar field. We define an analogous scalar field p that serves 
as the first field of the coordinate system, i.e., 

P=(^-C 2 )cos9. (12) 

The gradient Vp defines a vector field that points in the direction of the magnetic field (and 
hence points in the direction of the outflow). Next we construct the perpendicular vector 
field Vg, where the second scalar field q provides the second coordinate and is given by 

q={Pe + 2/0 V2 sin^ (13) 

The pair (p, q) thus represents a set of perpendicular coordinates in the poloidal plane, 
and these can be used instead of the original spherical coordinates (£, 9) or the cartesian 
coordinates (x,z). In this version of the problem, the field is axisymmetric about the i 
axis, so we can use the usual azimuthal coordinate as the third scalar field of the set 
(the magnetic field is poloidal, with no toroidal component). Note that both p and q are 
dimensionless; one can reinsert factors of the planetary radius Rp (as necessary) to convert 
back to physical units. Notice also that in the limit /3 — > 0, one recovers the coordinates for 
a dipole (Radoski 1967). 

The set of covariant basis vectors arises from the gradients of the scalar fields that 
define the coordinates. If we express these basis vectors in terms of the original spherical 
coordinates (£, 9, (/>), the basis takes the form 

e p =(/3 + 2r 3 ) cos 9 f - (/3 - T 3 ) sin 9 9 , (14) 

ig = (p + 2r 3 r 1/2 [(/? - r 3 ) ^ e ? + (p + 2c 3 ) cos e § , ] , (15) 

and 

^ = T —J- (16) 
v £ sin 

The quantities (f, 9, 0) are the usual unit vectors for spherical coordinates. However, one 
should keep in mind that the are basis vectors (not unit vectors), so that their length is 
not, in general, equal to unity (for further discussion, e.g., see Weinreich 1998). In general, 
the corresponding scale factors are given by the relation 

hi = _1 , (17) 
so that the scale factors for this coordinate system can be written in the form 

h p = [(P + 2r 3 ) 2 cos 2 9 + (P- T 3 ) 2 sin 2 9] ~ 1/2 , (18) 
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h q = (P + 2r 3 ) 1/2 [(/3 + 2r 3 ) 2 cos 2 + (/3 - r 3 ) 2 sin 2 0] _1/2 , (19) 

and 

= C sin . (20) 
The general form of the divergence operator is thus given by 



v.v 1 



hphqh^ 



d d 
g£ {hqh+Vp) + — (hph^Vg) 



4 sm 9 a<p 



In this problem, the fields are axisymmetric so that the derivatives vanish. Further, for 
flow along field lines, the vector fields (e.g., the velocity field) have only one component and 
depend on only one coordinate, so that the divergence operator collapses to the form 

^ hphghfj, dp ^ 9 * p ^ h p dp ~^ h v h q h^ dp ^ 9 ^ ^ 

For convenience, we define the following ancillary functions: 

f = /3 + 2C 3 , g = /3-C 3 , and H = f cos 2 9 + g 2 sin 2 9 . (23) 
In terms of these functions, we note that {h q h^) = qH~ 1 / 2 and that |B| 2 = HB 2 P . 

This specification of the divergence operator is implicit. One could invert equations 
f fT2|) and f lT3|) . and then write the spherical coordinates (£,9), the scale factors (h p , h q , h^), 
and the ancillary functions (/, g, H) as functions of the new coordinates (p, q). However, the 
definitions of (p,q) are cubic functions of (£, 9), so that the solutions of the cubic inversion 
are complicated and unwieldy (but still can be written down analytically). For clarity, we 
leave this construction in implicit form. 



3.2. Fractional Active Area of the Planetary Surface 

With the field configuration of equation (|lip . some fraction F^p of the magnetic field 
lines that originate on the planetary surface are open, whereas some field lines curve back onto 
the planet. Using this simplified model for the magnetic field, we can determine the fraction 
Fap- In this context, the open field lines that continue to "infinity" in this reduced problem 
join onto the magnetic field lines of the star and stellar wind (although the subsequent 
curvature of the field lines is not captured by this model). One can define a benchmark mass 
outflow rate to be that obtained for spherical flow over the entire planetary surface. We note 
that the fraction Fap of the surface that supports open field lines is related to — but is not 
equivalent to — the fraction of this benchmark outflow rate that the planet produces (see 
Section 4.4). 
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The first step is to solve for the magnetic field lines. For this field configuration, this 
construction can be done analytically. Assume that a field line begins at coordinates (£ = 
l,#o) corresponding to the planetary surface. The field line can then be represented as a 
curve in the plane such that 



sm9 J ^2 + /3£ 3 

Note that this expression is equivalent to the statement q = sin 9 (2 + (3) 1/2 = constant] this 
result must hold since the coordinate q was constructed so that Vg is perpendicular to the 
field lines. Using this solution, we can solve for the angular coordinate Ox for the field line 
that goes through the X-point. In this case, the X-point radius £x is given by the condition 

= 1 • (25) 

Note that this result is exact for the reduced field configuration considered here. The critical 
angle 9 X is given by 

sin 2 0x = 3/3 1 / 3 /(2 + /3) , (26) 

where this expression is valid for [3 < 1. The field lines that originate at small angles, 
< 9 < Ox-, with 9x defined above, are those that are open (reach spatial infinity) in this 
reduced problem. In the full problem, these open field lines join onto the background field 
lines of the star and stellar wind. The resulting fraction F^p of the surface area of the planet 
that can support an outflow is then given by 



F AP = 1 



3/3 1 / 3 
2 + f3 



1/2 



(27) 



Note that for /3 > 1, the background (stellar) field dominates that of the planet, field lines 
originating from all planetary latitudes are open and hence F^p = 1. In the opposite limit 
of small /3c 1, the fraction Fap ~ 0.75/3 1/3 ; since (3 ~ 10 3 for typical cases, only about 10 
percent of the planet surface can support outflow. 

With this field geometry, the field lines become asymptotically straight, so that the 
outflow has a cylindrical form in the limit £ — > oo. The outer surface of the resulting outflow 
cavity is delineated by the field line that passes through the X-point. If we use equation ( 1241) 
for the critical streamline, and then take the limit £ — > oo, we find that the radius of 
this cylinder is given by 

Woo = VSr x = V3(3- 1/3 R P = V3(B P /B*) 1/3 (aR P /R*) . (28) 



This cylindrical flow represents the outer limit of the inner problem, the regime where the 
outflow is launched from the planetary surface. This flow also represents the inner limit of 
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the outer problem, the regime where the flow follows the stellar field lines. Note that for 
typical parameters w 10 Rp ~ R*. 

Since the magnetic field strength vanishes along the critical streamline, the magnetic 
pressure will not be strong enough to dominate the ram pressure near the X-point. We can 
estimate the fraction of streamlines that are affected by this issue: Consider the surface 
where £ = £x = /3 -1 ^ 3 . For a given angle 6 on this surface, the magnetic field strength (from 
equation [IT] ) is given by B = 3Bp(3 cos 8. For the flow to remain magnetically controlled, 
the field strength must be larger than a critical value Bq determined by the ram pressure 
of the outflow (B c = 8npv 2 ), where we expect Be ~ 0.001 G for typical cases (see equation 
[5]). This condition requires cos 6 > Bc/(3Bp(3), which in turn restricts the coordinate q 
that specifies the streamlines to the range defined by 



Since we expect Be ~ /3Bp, the range of streamlines that are unaffected by the X-point issue 
is roughly given by q < qxV8/3. In other words, for ~94 percent of the streamlines, the flow 
remains magnetically controlled over the entire space. For the remaining ~6 percent, the 
flow is magnetically controlled except for a small region (of size (AL) ~ surrounding 
the X-point. Note that outflow can still take place along these streamlines, but that the flow 
pattern will be slightly different than that given by the unperturbed magnetic field geometry. 
For the remainder of this paper, for simplicity, we consider the flow to take place over the full 
range of open streamlines q < qx, while noting that this approximation may overestimate 
the outflow rates by a few percent. 



Given the specification of the magnetic field structure and hence the flow geometry 
(Sections [2] and [3]), we now consider the launch of a wind or outflow from the surface of the 
planet. After writing down the full set of equations of motion, we consider a reduced version 
of the problem where the flow is taken to be isothermal. For this case, the solutions for the 
dimensionless fluid fields can be found analytically, including the required conditions for the 
flow to pass smoothly through the sonic transition. In order to complete the solution, we 
must then specify the values for the physical parameters, i.e., the density p\ at the inner 
boundary and the sound speed a s (taken to be constant). 




(29) 



4. Outflows from Planetary Surfaces 
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4.1. Formulation of the Wind/Outflow Problem 



The equations of motion for this problem include the continuity equation, 



dp 
dt 



+ V-(pu) = 0, 



(30) 



the force equation, 




W - -VP + 

P 



Anp 



1 



(V x B) x B 



(31) 



the energy equation 




,(§ + ..v*)~ 

n for the magnetic fielc 



PV • u + r - A , 



(32) 



dB 

~dt 



+ V x (B x u) 



V x (r/ rcs V x B) 



(33) 



where ?7 res is the resistivity. This paper considers the magnetic fields to be fixed, and strong 
enough to not be changed by the outflow. In a full treatment, however, the back reaction of 
the outflow on the magnetic field should be taken into account. 

In this paper, we consider the gravitational potential ^ to be that of the planet, which 
is taken to be spherical with mass Mp and radius Rp. Since the planet orbits the star, the 
full potential has an additional contribution from the rotating frame of reference. The order 
of this correction term is 0(M*r 3 /Mpa 3 ), so that it has has size ~ 10~ 3 near the planet 
surface and size ~ 0.03 near the sonic surface. As a result, this term does not greatly affect 
the launch of the outflow and is not included here. To consistent order, we also ignore the 
tidal forces from the stellar gravitational field (note that the sonic surface is well inside the 
Roche radius). We thus work in the regime where the gravitational force is dominated by 
that of the planet. 

In the energy equation ( I3"2"j) . £ is the specific energy of the fluid, T is the heating rate 
(per unit volume), and A is the cooling rate. The heating is primarily due to UV flux from 
the central star (e.g., see MCM and references therein) so that T can be written in the form 



where ^ a b s is the fraction of the UV energy that is deposited by heat, er uv is the cross section 
for UV photons, F uv is the unattenuated UV flux of the star at the location of the planet, 
and r is the optical depth (from the star to the point where the heating term is evaluated). 



r = r] ahs F uv e T cr uv n , 



(34) 



- 16 - 



The cooling process is primarily due to Lyman-a radiation that is emitted by Hydrogen 
atoms as they are excited via collisions (Black 1981). To leading order, the cooling term can 
be represented by a function of the form 

A « Cn n + e- Tc/T , (35) 

where C = 7.5 x 10~ 19 , Tq = 118348 K; the resulting cooling rate has units of erg cm -3 
s _1 . Note that this particular form is only valid for temperatures T < 12,000 K; at higher 
temperatures the gas cools even more efficiently (Spitzer 1978). In any case, the cooling 
rate is large enough that the gas temperature never increases beyond an effective maximum 
T-max. ~ 10 4 K. In order for Lyman-a radiation to act as the primary cooling mechanism, it 
must dominate over other processes and the radiation must be able to escape; both of these 
conditions are met, as shown in the Appendices of MCM. 

Since the heating and cooling rates depend on the state of ionization, through the 
number densities of the neutrals n and ions n + , we also need the equation of ionization 
balance: 

+ V ■ (un+) = -jj^j 6 °"uv^o - oi R n + , (36) 

where the recombination coefficient o.r ~ 2.7xlO _13 (T/10 4 K) -0 9 cm 3 s _1 (Storey & Hummer 
1995). For the remaining parameters, we take cr uv = 2 x 10~ 18 cm 2 and (hu) = 2.2 x 10 11 
erg (e.g., Spitzer 1978). 

Note that the planet is only heated on the side facing the star, whereas cooling takes 
place over the entire surface of the planet. This work implicitly assumes that the UV heating 
is distributed uniformly throughout the upper atmosphere of the planet. This approximation, 
in turn, is valid when zonal winds are strong enough to provide the required redistribution. 
Although the issue is not settled, current models suggest that strong winds are present, so 
that uniform heating is a reasonable approximation (for further discussion, e.g., see Batygin 
& Stevenson 2009, Langton & Laughlin 2008, and references therein). 

As formulated here, equations (1301-1361) make up a compete set that can be solved for the 
fluid fields. However, even if we use a simplified magnetic field geometry, with the coordinates 
constructed in Section [3j the problem remains intrinsically three dimensional: Although the 
flow is axisymmetric and follows the coordinate p defined by equation (1121) , the heating comes 
from the central star which lies off to one side. As a result, the incoming photons do not 
follow the coordinates, and a full solution for the heating/ cooling of the outflow requires one 
to solve a three-dimensional radiative transfer problem. Before embarking on that task, it is 
useful to have solutions for an approximate treatment. Toward this end, a simplified version 
of the problem is formulated in the next section. 
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4.2. Reduced Equations of Motion 

In this section, we consider steady-state solutions and assume that the magnetic field 
structure due to the planet (and the star) are strong enough to dominate the flow. As a 
result, in this regime, the magnetic field is fixed and current-free. The continuity, force, and 
induction equations thus reduce to the forms 

V • (pa) = , u ■ Vu + W + - VP = , and B = Tpu , (37) 

P 

where the parameter T is constant along streamlines (e.g., Shu et al. 1994, Cai et al. 2008). 

The velocity vector u follows the magnetic field lines, which follow the coordinate p 
in the system constructed in Section [3J In other words, the flow velocity has only one 
component, which points in the direction of the magnetic field p = h p e p (by construction). 
Next we assume that the flow is isothermal with constant sound speed a s and define the 
following dimensionless quantities 

lul or \I/ 

u = u(p) = — , a = —, C = and ^ = -» • (38) 

a s pi Hp aj 



Here, Rp is the radius of the planet and pi is the density at the inner boundary £ = 1. The 
continuity equation thus takes the form 

du da au d , . , . 

% + M ^ = -V^ (W ' (39) 

and the force equation becomes 

du ^ 1 da dip dip <9£ 
dp a dp dp <9£ dp 

These equations can be integrated immediately to obtain the solutions 

auhqh^ = auqH" 1 ^ 2 = X , (41) 

and 

-u 2 + log a + ip = e . (42) 

Although the potential ip, in general, contains additional contributions (e.g., tidal forces), 
we specialize to the case that includes only the gravitational potential of the planet so that 
ip = — &/£, where b = GMp / \a 2 s Rp). Note that the quantity h q h<j, = qH^ 1 ^ 2 is proportional 
to the inverse of the magnetic field strength (consistent with the third part of equation [37]). 
The parameters A and e are constant along streamlines, but are not, in general, the same 
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for all streamlines (they are functions of q). In order for the flow to pass smoothly through 
the sonic point, only particular values of the constant A are allowed. This constraint is 
considered in the following section. 

The boundary conditions at the planetary surface take the form 

£ = 1 , a = 1 , and u — U\ — XH^ 2 jq . (43) 

Since A is determined by the conditions at the sonic point, u\ is specified. In addition, the 
remaining parameter e is determined by evaluating the force equation at the inner boundary, 
i.e., 

e=\u\-b = ¥^-b. (44) 



The outflow starts with subsonic speeds so that ui < 1 (below we find that u\ ~ X/q ~ 
0.01), whereas typical planet properties imply that b ~ 10. As a result, one can use the 
approximation e ~ — b with good accuracy. 



4.3. Sonic Point Conditions 

Critical points in the flow arise when the fluid speed is equal to the transport speed. 
In general, magnetic media support three types of MHD waves and hence allow for three 
types of critical points (e.g., Shu 1992). In this case, however, the flow is confined to follow 
the magnetic field lines, so that only one possible critical point arises, in this case where the 
flow speed equals the sound speed. For the equations of motion f l4T|) and ( j4"2"j) . the required 
matching conditions at the sonic point take the form 



u 



2 



Id <9# d£ b d£ 

and j^^ {hM = ^ P = edp- (45) 

We must thus evaluate the geometrical factor Q defined by 

s =wM M - (46) 

Note that in spherical coordinates, this factor would have the usual form 2/£. In terms of 
the ancillary functions from equation ( 123]) . the partial derivatives can be expressed as 

d£ = fcos6 89_ = gcosO 

dp H an dp iH ' 1 } 
and the geometrical factor Q takes the form 

S = ^ [V 2 cos 2 9 - g(g + 2f) sin 2 9] . (48) 
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The matching condition at the sonic point can then be written in the form 

3 



[2f 2 cos 2 6-g(g + 2f)sm 2 6] . (49) 



One can eliminate the explicit angular dependence from this expression using the result 
sin 2 6 = q 2 /£, 2 f, and noting that q is a constant along the direction of the flow. The right 
hand side of the above equation (for given q) thus becomes a function of £ only, 

b = 2f 2 - (g 2 /f + 2g + 2f) q 2 ^ 2 

3 Pe + {9 2 -f 2 )q 2 ' 1 ] 

For given planetary properties (set by the value of b), magnetic field strength ratio (set by 
f3), and starting angle 9 of the streamline (set by q), equation (1501 provides an algebraic 
expression that can be solved for the value of £ = £ s at the sonic point. With £ s specified, 
the angle #(£ s ) is also specified, and hence the value of p = £gcos8 is determined. Finally, 
the value of the parameter A that allows for smooth flow through the sonic point is given by 



A = X(q) = qH; 1 / 2 



exp 



X 2 Hi b 1 
~W + Is 2 



(51) 



where the subscript 's' (T) implies that the quantity is evaluated at the sonic point (inner 
boundary). Equation ( l5Tj) provides an implicit solution for the parameter A. However, the 
A 2 term on the right hand side of equation (I5TI) is extremely small (it is equal to uf/2 <C 1) 
and can be ignored to leading order; doing so results in a direct expression for the parameter 
A (after equation [50]) has been solved to find the value of £ s ). We also note that one can 
define an alternate parameter A = X/q, which is useful because it is easier to find solutions 
for A (e.g., when q — > 0). 

Figure |2] shows the sonic surface for a planet / star system with typical properties. Here, 
the parameter /3 = 10~ 3 , which holds when the surface fields on the planet and star are 
equal, and the semimajor axis a = 10 i?*; the parameter b = 10, which holds for planets 
with mass Mp = 0.75M/, radius Rp = 10 10 cm, and sound speed a s = 10 km/s. The sonic 
transition occurs at nearly constant dimensionless radius (£ ~ 3.3), as delineated by the 
heavy curve. The light solid curves show the underlying coordinate system, i.e., the lines of 
constant p and q. Only the streamlines (constant q) corresponding to open magnetic field 
lines are shown. The last coordinate line shown is the one passing through the X-point. For 
the perpendicular coordinate, the value p = corresponds to the surface passing through 
the X-point, and that surface is a sphere. Negative values of p correspond to surfaces that 
are more highly curved (compared to a sphere), whereas positive values of p produce flatter 
surfaces. Note that the sonic surface lies well within the surface passing through the X-point. 
As a result, the flow diverges significantly faster than that of a spherically symmetric wind. 
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Fig. 2. — Sonic transition surface for planet with dimensionless gravitational potential depth 
b = GMp/(a 2 s Rp) = 10 and background magnetic field strength parameter /3 = 10 -3 . The 
heavy solid curve near £ = 3.3 shows the location where the flow passes through the sonic 
point. The light solid curves show the underlying (p, q) coordinate system. The cartesian 
coordinates x and z are given in units of the planetary radius. 
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Fig. 3. — Values of the parameter X/q that allow the outflow to pass smoothly through 
the sonic point. The parameters A are shown here as a function of sin^o = q(2 + /3) -1 ^ 2 , 
which determines the starting point of the streamlines. All cases use b = 10 (where b = 
GM P /(a 2 s Rp) sets the depth of the gravitational potential of the planet). Curves are shown 
for a range of the parameter that sets the relative strength of the stellar and planetary 
magnetic fields: (3 = 0.008 (bottom curve), 0.004, 0.002, 0.001, and 0.0005 (top curve). For 
smaller values of (3, the curves converge toward a well-defined locus (although the range in 
sin 6*o shrinks). 



-22 - 



The right hand side of equation (1501 is not, in general, a monotonic function of £; it 
reaches a maximum and then decreases in the limit £ — > oo. If the maximum value is too 
small, then equation ( 150]) has no solutions, and the outflow will not pass through the sonic 
point within the context of this simplified model. One can show that the right hand side 
of equation (150|) is a increasing function of q for sufficiently small £ (in the regime where 
matching occurs). As a result, the minimum value occurs for the smallest value of q, i.e., 
along the pole where q — 0. For this streamline, the matching condition is the most difficult 
to meet; the sonic point is thus given by solutions to the cubic polynomial b = 6//£ 2 , or, 
equivalently, 

/3£ 2 + 2/£ = 6/6. (52) 

This equation has no real solutions if either parameter ft or b is too large; solutions require 
ftb 3 < 8, or equivalently ft < /3 max = 8/b 3 . Since b ~ 10, the value of /3 max ~ 0.01. Only 
smaller values of ft < /3 max allow for smooth outflow solutions. 

For solutions that pass through the sonic point, so that A is specified, Figure |3] presents 
the resulting values of X/q over the range of allowed streamlines. These results are shown 
for a range of the magnetic field strength parameter ft. For the planetary properties used 
here, where b = 10, the maximum allowed value of ft — 0.008 is shown by the lowest curve 
in Figure |3j Larger values of ft do not allow for a sonic transition over the full range of open 
field lines. The resulting X/q for smaller values of ft converge toward the result for ft = 0. In 
addition, the range of allowed streamlines, given here by the range of allowed starting angles 
9 , decreases with decreasing ft (see equation [2B]). Notice also that X/q is relatively slowly 
varying over the range of allowed streamlines. For the largest possible ft, X/q varies by a 
factor of ~ 2 from the pole to the angle of the last open field line; the results for smaller ft 
show much less variation. 

Note that the results derived above make sense in the limits: For flow near the planetary 
surface, where the magnetic field is determined by the dipole of the planet, the sonic point 
condition of equation (|50|) reduces to the form b = 3£ (along the pole). This solution results 
from the continuity equation d^(l; 3 pu) = 0, where this form for the divergence operator is 
expected for flow that follows a dipole field. In the opposite limit where £ — > oo, the field 
lines and the streamlines become asymptotically straight and point in the z direction. In 
this limit, we recover the results for one-dimensional flow, where the continuity equation has 
the form d z (pu) = 0. The solutions that smoothly pass through the sonic point are those 
where the magnetic field of the planet dominates, so that the sonic point £ s ~ b/3. Using 
this result, we can define a fiducial value for the parameter X/q, i.e., 

b 3 

(A/g) = — exp [5/2 - b] « 0.2256 . (53) 
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For expected planetary properties with b « 10, we find (A/g)o ~ 0.01. Notice also that X/q 
sets the scale for the mass outflow rate and that this quantity decreases exponentially with 
increasing depth of the gravitational potential well of the planet (set by b). 

For solutions that smoothly pass through the sonic transition, we can also find the 
asymptotic speed Woo, i.e., the value realized in the limit £ — > oo. In this regime, the reduced 
equations of motion (I4ip and (142 p are simplified further to take the forms 

ua = Xj3/q and -u 2 + \oga = e, (54) 

where both A and e are known functions of the parameter q that specifies the streamline. 
The asymptotic speed u = is thus given by the expression 

u^-logi& = 2e + 21og(g//3A). ( 55 ) 

Note that equation (I55p has two roots whenever real solutions exist, which requires the right 
hand side to be greater than unity. One root corresponds to u 2 ^ < 1, whereas the (physical) 
root of interest corresponds to u 2 ^ > 1. 

Figure H] shows the dimensionless velocity and density profiles for outflows characterized 
by different values of the constant A. These profiles are shown for flow from the poles, so 
that q = and the flow direction is given by r = z. In the top panel, the central curve 
shows the velocity profile for the critical value of X/q that allows the flow to pass smoothly 
through the sonic point. If the value of X/q is too small (bottom curve), then the outflow 
speed never reaches the sound speed. Instead, the flow velocity reaches a maximum at the 
radius of the sonic point and then decreases for larger £. On the other hand, if the value of 
X/q is too large (top curve), the outflow cannot reach the sonic point in a smooth manner. 
This plot is thus analogous to that found for the well-known Parker wind (see Figure 1 of 
Parker 1965) and for Bondi-Hoyle accretion (e.g., see Shu 1992 for further discussion). 

The bottom panel of Figure H] shows the corresponding density profiles. Here the solid 
curve shows the result for the critical value of X/q that allows for a smooth sonic transition. 
The supercritical and subcritical cases are shown as the dashed curves. Note that the density 
profile is extremely steep near the planetary surface, and then levels out at large radii £. 
This behavior is a reflection of the divergence operator, which follows the magnetic field 
lines, which in turn spread rapidly near the planet and become straight at large distances 
£. Notice also that the density for the subcritical case is larger than that of the critical case 
for £ 3> 1. Although the mass outflow rate is smaller for the subcritical case (which lowers 
the density), the asymptotic speed (see the top panel) is much smaller and this latter effect 
increases the density and produces the observed behavior. 
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Fig. 4. — Top panel shows the velocity profiles for outflows from planetary surfaces for flow 
with q = (from the planetary pole). The three curves show the velocity as a function of 
radius £ for the critical value of X/q (central curve), for a subcritical case (bottom curve) 
and a supercritical case (top curve). These velocity profiles were calculated using b = 10 and 
f3 = 10~ 3 . Bottom panel shows the corresponding density profiles. The critical case is shown 
as the solid curve, whereas the dashed curves show the subcritical case (which continues to 
large £) and the supercritical case (which is confined to small £). 
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For completeness, we can find the values of the parameter X/q for field configurations 
where the matching condition (equation |49j ) at the sonic point cannot be met. In this case, 
the largest value of A, and hence the largest outflow rate, occurs when the outflow velocity 
approaches the sound speed in the limit £ — > oo. In this case, we can combine equations 
( 14"T|) - (jHJ), while setting u — > 1, and then solve for the value of the parameter X/q: 

/T^xp [-(6+1/2)] . (56) 



X/q = f3 1 exp 



-(6 + 1/2) + 



A 2 ill 

2q 2 



4.4. Dimensionless Mass Outflow Rate 

In this geometry, the continuity equation reduces to the form auqH' 1 / 2 = X = constant 
(see equation |4T)). Because the quantity qH~ 1 / 2 does not have units of area, the constant 
A is not the mass outflow rate. As a result, we need to find the relationship between the 
outflow rate and the quantities that appear in the equations of motion. 

In the limit of large £ — > oo, the dimensionless mass outflow rate is given by the 
integral 

riioo = 2 / 2-Kwdw (cm)oo , (57) 



Jo 

where w = w/Rp and where the dimensionless outer radius of the outflow can be 
determined from equation (|28p . The leading factor of two arises because the wind flows 
from both the northern and southern hemispheres of the planet. From equation (|4ip we 
find that (au)^ = H l l 2 u\H 1 1//2 , where the subscripts '1' indicate that the quantities are 
evaluated at the inner boundary. In the limit £ — > oo, f — > {3, g — > /3, and H — > [3 2 . Using 
the streamline equation, the cylindrical radius w is related to the orthogonal coordinate q 
through the expression w 2 = q 2 / '/3. After changing variables, the dimensionless outflow rate 
takes the form 

fix 

rhoo = 471 U\qdqH l , (58) 



o 



where Hi is a function of q and where q\ = 3/3 1 / 3 . In general, the starting speed U\ will not 
be the same for all streamlines, i.e., it will be a function of q. If we replace u\ in the integral 
by the appropriate mean value (ui), equation f )58|) can be evaluated to find 

9^/3(1 + 2/3) i 1/2 ' 



°° ~ 3(1 + 2(3) ' 



1 - 



(2 + /3Y 



(59) 



This expression is valid for (3 < 1. For larger values of (3, the stellar contribution dominates 
the magnetic field of the planet over its entire surface, and the dimensionless mass outflow 
rate reduces to the spherically symmetric form rh = Aitui. 
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The mass outflow rate can also be evaluated at the inner boundary £ = 1. The mass 
outflow rate must be the same in both limits £ — > 1 and £ — > oo, so this calculation provides 
a consistency check. In this case, the dimensionless outflow rate is given by the expression 



where the integral is taken over the entire planetary surface; however, the outflow velocity 
u is only nonzero over the fraction of the surface that supports the outflow. Notice also 
that the flow is not radial, so that the velocity u points in the p direction (rather than the 
f direction). Using the divergence theorem, the integral can be rewritten in the form 



where V is the volume of the planet and where u is the only non- vanishing component of the 
velocity (in the p direction). Here, the volume element dV = h p h q h^dpdqd(j). The integral 
with respect to <ft produces a factor of 27r because the system is axisymmetric; the integral 
over p can be evaluated directly and results in a surface term evaluated at the planetary 
surface £ = 1. The remaining expression becomes 



J o 

which has the same form as that in the limit £ — > oo (compare with equation [58]). Again 
we must consider flow from both the northern and southern hemispheres of the planet. 

Figure [5] shows the dimensionless mass outflow rate fa = dm/dt plotted as a function 
of the dimensionless depth b of the gravitational potential well of the planet. The results 
are shown for a series of values of the field strength parameter /3; the curves correspond 
to increasing values of j3 from bottom to top in the figure. Notice that the curves end at 
particular values of b, and that these values decrease with increasing 0. For larger values 
of b, only a fraction of the open streamlines allows the flow to pass smoothly through the 
sonic point. The mass outflow rates are thus diminished. In the discussion below we derive 
a scaling law (see equation [64]) that describes how the outflow rate m depends on the 
variables (b, over the parameter space represented in Figure El this scaling law holds to 
an accuracy of ~ 3%. 

Figure [6] provides another way to view the dependence of the dimensionless outflow rates 
on the underlying parameters. In this case, we plot the outflow rate m as a function of the 
field strength ratio /3. Each curve corresponds to a different value of the dimensionless depth 
b of the gravitational potential well of the planet. The outflow rates increase with /3 until a 




(60) 




(61) 




(62) 
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Fig. 5. — Dimensionless mass outflow rate versus dimensionless depth b of the gravitational 
potential well of the planet. The curves show the results for different values of the magnetic 
field strength parameter j3, from j3 = 0.0005 (bottom curve) to j3 — 0.008 (top curve). 
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Fig. 6. — Dimensionless mass outflow rate versus magnetic field strength parameter (3. The 
curves show results for different values of the dimensionless depth b of the gravitational 
potential well of the planet, where b = 12 (bottom), b = 10 (middle), and b = 8 (top). Local 
maxima occur for the largest values of /3 that allow the flow along all of the open streamlines 
to pass smoothly through the sonic point (see text). 
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well-defined maximum is reached; note that the parameters of this extremum are defined by 
solutions to equation fl52|) . At this point, further increases in (3 lead to straighter magnetic 
field lines (and streamlines), especially along the poles of the system, and the flow along 
those streamlines cannot pass smoothly through the sonic point. As f3 increases, the fraction 
of the open streamlines that allow for smooth flow decreases, but the fraction of streamlines 
that are open increases. In addition, the sonic point moves inward, for those streamlines 
where it exists, and this effect acts to increase the outflow rate rh. These competing effects 
thus lead to the non-monotonic behavior shown in Figure El 

The parameter is typically small; for example, ~ 10~ 3 when the field strengths are 
equal on the stellar and planetary surfaces and a = 10 R* ~ 0.05 AU. As a result, it is useful 
to find simplified results that are correct to leading order in 0. If we expand equation ([59]) . 
the leading order term becomes 

m = 3Tt(u 1 )p 1 / 3 + O(0 2/3 ) . (63) 

Note that in the case of a spherically symmetric flow, the dimensionless mass loss rate would 
have the form m sp h = Anui in these units. The fraction of the total possible outflow rate (with 
fixed Mi = (ui)) is thus F ~ 3/3 1//3 /4. From the previous section, the fiducial value of the 
parameter (A/g) ~ U\ ~ 6 3 exp[— b] (see equation [53]). We thus except the dimensionless 
mass outflow rate rh to obey the scaling relation 

rh w A m b 3 exp[-b] 1/3 , (64) 

where A m is a constant of order unity. Fitting to the results presented in Figures Eland El we 
find A m ps 4.8 ±0.13 (where the quoted uncertainties represent the standard deviation of A m 
for the parameter space depicted in the Figures). This scaling law works well as long as the 
flow along all of the open field lines can pass through the sonic point. For sufficiently large 
0, however, sonic transitions cannot take place in the polar direction, and the dependence 
of the outflow rate rh on the parameters (b, 0) becomes more complicated. 

In the limit of large 0, sonic transitions cannot take place along any of the directions. 
On the other hand, all of the field lines from the planet must match onto stellar field lines. 
In this case, the values of X/q are given by equation (1561) for all of the streamlines. The 
integral that defines the dimensionless outflow rate can be evaluated to obtain the result 

m = 6vr/r 2/3 exp [-(& + 1/2)] . (65) 

This expression is only valid for < 1. At larger stellar field strengths, all of the field 
lines originating on the planet surface are open, and the leading coefficient in equation (165]) 
becomes 27r(l + 2/0). In the extreme limit ^> 1, the dimensionless outflow rate thus 
reduces to the expression rh = 27rexp[— (6+ 1/2)], which is the form expected for flow in the 
z direction from a disk with the radius of the planet. 
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4.5. Estimating the Physical Constants 

The previous subsections specify the solutions for the dimensionless fluid fields, including 
the necessary conditions for passing smoothly through the sonic point and specification of the 
dimensionless mass outflow rate m. In this section, we complete the solution by estimating 
values for the physical parameters p\ and a s that determine the full mass outflow rate, where 
M = pia s Rpm. 

We first provide order of magnitude estimates: The planetary radius is given, with 
a typical value Rp ~ 10 10 cm. Since these planets are gas giants, they do not have solid 
surfaces. As a result, the planetary radii measured by transits are determined by the levels in 
the atmosphere that are opaque to optical photons, typically at pressures ~ 1 millibar (see, 
e.g., Charbonneau et al. 2007). We also note that the launching radius where £ = 1 could 
lie several scale heights above the planetary surface defined by transit measurements (e.g., 
at pressures ~ 1 nanobar). However, the scale height H « Rp/b, where b ~ 10, or larger. 
As a result, H -C Rp, and any departures of the launching radius from R P are expected to 
be small (~10 percent). The gas temperature is expected to approach the benchmark value 
T ~ 10 4 K, so that sound speed a s ~ 10 km/s. Finally, the density pi, or equivalently rii = 
pi/mp, can be estimated by using the fact that the wind is launched near the r = 1 surface. 
The optical depth r ~ a uv riiH, so the number density n\ ~ l/(Ha uv ) ~ b/(Rpcr uv ) ~ 10 9 
cm -3 . With these values, the scale for the mass outflow rate is M sca i e = mpnia s Rp ~ 10 11 
g/s (see also MCM). Since the dimensionless factor fn ~ 0.01 (see Figures |5] and |6l and 
equation [M]), typical outflow rates are expected to be M ~ 10 9 g/s. 

To obtain a better estimate for the quantities (a s ,rii), we need to consider heating 
and cooling of the gas (equations [22], [31], and [25]) including ionization (equation [3"o]). 
Unlike previous outflow studies, the geometry of the flow is determined by the magnetic 
field structure. For relatively large (3, the magnetic field of the star dominates, and the 
streamlines become primarily vertical. In the opposite limit of small (3, only the streamlines 
from the polar regions lead to outflow, and these streamlines are also oriented mostly in 
the z direction (see Figure [2]). In either case, the flow is directed in the polar directions, 
whereas stellar heating arrives from the equatorial direction. Because of this configuration, 
the optical depth of the incoming radiation is not directly tied to the mass outflow rate (which 
occurs for spherical flow). In addition, the stellar UV photons penetrate (from the side) into 
the atmospheric layers where the flow speeds are small. Leaving a full three-dimensional 
treatment of the heating/cooling problem for the future, we adopt here a simplified approach 
where flow velocities are neglected. As a further approximation, we assume that the required 
values of the sound speed a s and density ri\ are determined at the layer where the optical 
depth is unity, where n x = pi/m P is the total number density at the r = 1 surface. 
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In this limit, the energy equation (132]) reduces to the condition that heating and cooling 
are locally in balance, so that A = T. Similarly, the ionization equation ( 136]) reduces to the 
statement that the rate of ionization balances the recombination rate. These two equations 
can be combined to determine the ionization fraction as a function of temperature. With 
the definition X + = n + /rio, we can write 

X+(T) = ?± = expl-T c /T} « (2.5 x 10 5 ) T 4 a9 exp[-11.8348/T 4 ] . (66) 

With the ionization fraction specified, the heating equation determines the temperature and 
can be written in the form 

eXp[ - Tc/T] ITx; = Cn eXp[ ' r]= Ce ^ ' (67) 

where n = no + n + is the total number density and n\ = n{r = 1). The left hand side 
of equation flSTj) is a function of temperature, whereas the right hand side is a function of 
position only. The second equality specializes to the layer where r = 1. If we assume that 
the flow velocities are small, the optical depth is given by 

Tiicr uv R p i , kq\ 

where b is the dimensionless depth of the gravitational potential well (and note that b oc 1/T). 
Equations (IBTj) and (!68|) provide two equations for the two unknowns n x and T. With these 
parameters specified, we can then evaluate the scale M sca i e = mpriia s Rp for the mass outflow 
rate. The outflow rate itself is given by M = m(mpnia s Rp) , where the dimensionless outflow 
rate fn can be approximated using the scaling law from equation or the full calculation 
in that section. Note that in order to determine the value of b = GM P / a 2 s R P) we need the 
temperature T oc a 2 s . 

Figure [7] shows the number densities n\ and corresponding outflow rates M as a function 
of the UV flux Fjjy for three choices of planetary mass: Mp = 0.5, 0.75, and 1.0 Mj. Keep in 
mind that the dimensionless outflow rate m depends on the magnetic field strength parameter 
(3 (see Figure [6]); these curves are calculated using (3 = 0.001. In this approximation, the 
outflow rate increases somewhat more slowly with F\jy than the linear relation of equation 
([1]). For sufficiently large flux levels (not shown), the temperature approaches its effective 
maximum value (just above 10 4 K, see Spitzer 1978); in this regime, the outflow rates would 
approach a constant value, but the assumptions of this section break down. Notice also that 
the number density (dashed curves) is a slowly varying function of both the UV flux and the 
planet mass, and that rii ~ 10 9 cm -3 as expected. 
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Fig. 7. — Mass outflow rates as a function of UV flux from the star. The outflow rate dM/dt 
= M is given in g/s, and the UV flux is given in erg cm" 2 s _1 . Solid curves show the outflow 
rates for three values of the planet mass, from Mp = 0.5 Mj (top curve) to Mp = 1.0 Mj 
(bottom curve). Dashed curves show the number density n\ (in cm" 3 ) at the r = 1 surface 
for the same cases. The magnetic field strength parameter f3 = 0.001 for all cases shown. 
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5. Observational Signatures 



Observations used to infer the presence of planetary outflows show that the transit depth 
is greater at UV wavelengths than in the optical (Vidal-Majar et al. 2003, 2004; Desert et 
al. 2008, Sing et al. 2008, Lecavelier des Etangs et al. 2008, Linsky et al. 2010). As a result, 
we need to determine the optical depth of the outflow to UV radiation. To start, we define 
the dimensionless column density iV c according to 



where a is the dimensionless density obtained from the flow solution (see Section H]). The 
column density N c is thus defined for a given path and the variable s is the dimensionless 
distance along the path. 

Note that equation (1691) determines the total column density, whereas the optical depth 
depends on the chemical species that absorbs the UV radiation. Most observations to date 
are carried out for Lyman-a photons, so that the column density of interest is that of 
neutral hydrogen. As a result, ionization modeling is necessary to determine the optical 
depths and hence the observational signatures. This present treatment assumes isothermal 
flow. To consistent order of approximation, the neutral column density is reduced from that 
of equation ( 1691) by a factor F = (1 + X + ) _1 (see equations [67] and [68]). For the benchmark 
temperature T = 10 4 K, X + « 1.8, so that the neutral hydrogen column density is reduced 
to a fraction F pa 0.36 of the total column density. However, the ionization fraction is a 
sensitive function of temperature. For T = 8000 K, X + pa 0.08, and the factor F f» 0.93. 

For the sake of definiteness, we evaluate the column density on paths that are parallel 
to the star-planet direction (y = in the coordinate system used here; see Figure [Q with a 
constant vertical height z. Each point along the path corresponds to a different streamline, 
and hence a different value of the constant X/q. Figure [8] shows the dimensionless column 
density N c through the outflow, as a function of z, for several typical cases. The solid curves 
show the effect of varying the depth of the gravitational potential well of the planet. The 
three curves correspond to b = 8, 10, and 12, where the magnetic field strength parameter 
(3 = 0.001, which corresponds to equal surface fields on the star and planet (for a = 10 R*). 
The dashed curve shows the effect of doubling the field strength parameter (3 (for b = 10). 

The optical depth at UV wavelengths is given by r = Friia nv R P N c . For typical cases, 
the quantity riia UY Rp ~ 20. The base density nj = pi/m P is a slowly increasing function 
of the incident UV flux Fuv ( see Figure [7J, but the correction factor F due to ionization 
decreases with temperature and hence decreases with increasing Fuv- The range in z for 
which the flow is optically thick is thus relatively small, ranging from z/Rp = 1.35 to 1.77 
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Fig. 8. — Dimensionless column density through the outflow as a function of vertical coordi- 
nate z (in units of the planet radius Rp). Solid curves show the column density for magnetic 
field parameter ft = 0.001 with three values of the dimensionless depth b of the planetary grav- 
itational potential well: b = 8 (top), 10 (middle), and 12 (bottom). The dashed curve shows 
the b = 10 case with larger field strength parameter /3 = 0.002. Optical depth is obtained 
from the dimensionless column density through the relation r = Fnia^RpNc ~ 20iV c . 
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for the cases shown with b = 8 to 12. With these values, the area for which the planet is 
apparently optically thick will be 2 - 3 times larger at UV wavelengths than in the optical 
bands. The transits for UV observations are thus predicted to be 2 - 3 times deeper. 

The mass outflow rates are expected to vary substantially from planet to planet. In 
order to illustrate this trend, in Figure [9] we plot estimates for the mass outflow rates for the 
collection of extrasolar planets that are observed in transit (data from Schneider 2010). The 
planetary masses, radii, and semimajor axes are observed. In order to estimate the outflow 
rates, we must also specify the UV flux F uv and the magnetic field strength parameter (3. 
With these quantities determined, the results of Section @] can be used to calculate the mass 
outflow rate M. The UV fluxes from main-sequence stars vary with spectral type and other 
stellar parameters (e.g., rotation rates) in a complicated manner (see Lecavelier des Etangs 
2007, Lammer et al. 2003, and references therein). To obtain the results shown in Figure El 
we estimate the UV flux using a simple scaling law that is intermediate between the scaling 
laws advocated by the aforementioned authors. To specify the magnetic flux parameter (3, we 
assume that the stellar and planetary surface fields are comparable, so that to leading order 
oc a~ 3 . Finally, we take into account the fact that radius of the r = 1 surface (where the 
outflow is launched and the parameter b = GM P /a 2 s R is evaluated) lies above the planetary 
surface at Rp. Using a simple hydrostatic model for the lower layers of the outflow region, 
we estimate that R\ ~ 1.2Rp, and this correction is used here. 

With the specifications described above, we obtain the mass outflow rates shown in 
Figure M as a function of planet mass. Several trends are clear: First, for a given planet 
mass, a wide range of outflow rates are possible. Both the magnetic field parameter (3 and 
the UV flux Fuv vary with semimajor axis a, so that closer planets are predicted to have 
stronger outflow rates. In addition to this dependence on a, however, both variables (/3, F uv ) 
are expected to display substantial variation from system to system. Next we note that range 
of outflow rates shown in Figure M spans seven orders of magnitude. Those planets with the 
largest outflow rates should thus have observable transit signatures, whereas those on the 
other end of the range will show little or no effect. Finally, the outflow rates are predicted to 
be a steeply decreasing function of planetary mass Mp. In the isothermal limit, we obtain 
the nearly exponential decrease shown in Figure M (see equation [M])- If the gas can cool 
substantially before it passes through the sonic transition, this decrease could be less steep, 
but the overall trend remains. 

Given that two observed planets have signatures of planetary outflows, it useful to see 
how they compare with the results shown in Figure IHJ The planet HD209458b (Vidal-Majar 
et al.2003) has mass mp = 0.685 mj, and the figure shows that the expected outflow rates 
are M ~ 10 10 g/s. Even with the reduction due to magnetic effects, such estimated outflow 
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Fig. 9. — Rough estimates for mass outflow rates for the sample of extrasolar planets observed 
in transit. Mass outflow rates dM/dt = M are given in g/s; planetary masses are given in 
Jovian masses Mj. 
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rates are large enough to account for the observations of HD209458b (although these values 
are near the low end of the inferred range). However, for the planet HD189733b (Lecavelier 
des Etangs et al. 2010), the mass is larger, m P = 1.15 mj, and the predicted outflow rate 
from this theory is much lower. The values shown in Figure M correspond to M ~ 10 8 g/s 
with significant scatter. The observational papers for HD189733b advocate a large UV flux 
(up to 40 times the solar value), which increases the theoretical outflow rate to M ~ 10 9 
g/s for our choice of field strength parameter /3 = 10~ 3 . To account for the inferred outflow 
rate of M ~ 10 10 g/s, one could invoke a larger value of j3 for this system. Nonetheless, 
some tension remains between the inferred outflow rate and that expected theoretically from 
magnetically controlled models. More observations, of this system and others with a range 
of masses, are need to sort out this comparison. 



6. Conclusion 

6.1. Summary of Results 

This paper has begun a theoretical study of outflows from the surfaces of Hot Jupiters 
in the regime where the flow is controlled by magnetic fields. In this case, the magnetic field 
structure determines the flow geometry (whereas the field configurations are determined 
by independent dynamo processes within the planet, and by the background contribution 
from the star). With the magnetic field structure specified, the dimensionless version of the 
outflow problem can be solved semi-analytically; this paper carries out this calculation in 
the isothermal limit, including the requirement that the flow must pass smoothly through 
the sonic transition. The determination of the physical constants represents the final piece 
of the calculation. The specific results of this paper can be summarized as follows: 

[*] This paper considers a reduced description of the magnetic field structure that in- 
cludes the dipole field of the planet and a constant background contribution. The reduced 
field structure near the planet is modeled in Section [3l and provides a description of the 
flow geometry in the region where the outflow is launched (see Figure [2]). In this section 
we construct the corresponding coordinate system and differential operators for this flow 
geometry. These results are used here to study the launch of outflows for Hot Jupiters, but 
can also be used in a variety of other applications. 

[*] A fraction Fap of the planetary surface supports magnetic field lines that are open 
with respect to the planet. This area fraction provides a constraint on the mass flow that can 
fully leave the planet. The remaining fraction, 1 - F^p, defines the region where material 
can leave the planetary surface but is nonetheless confined to the immediate vicinity of the 
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planet; this material produces an exosphere surrounding the planet. In the magnetically 
controlled regime, the outflow rate is thus lower than estimates obtained by assuming spher- 
ical symmetry. The fraction F^p is calculated analytically for the reduced field configuration 
near the planet (Section [3]). 

[*] Along each streamline, the dimensionless energy e and the flow momentum parameter 
A are constant, but vary across streamlines. These parameters are specified by requiring 
that the flow pass smoothly through the sonic point, and by the boundary conditions at 
the planetary surface. The resulting specification of parameters can be found analytically 
(see equations [50] and [51] . and Figure [3]). The outflow rate is determined by integrating 
over the surface area of the planet where the outflow is active, where the angle of the flow 
direction (which is not radial) must be included (see equations [33 T ne resulting mass 

outflow rates are well-defined functions of the depth b of the planetary potential well and the 
magnetic field strength parameter (3 (see Figures [5J and [6]) . Over the regime of parameter 
space where flow along all of the open field lines can pass smoothly through the sonic point, 
the scaling law of equation ( 164|) provides an accurate determination of the dimensionless 
outflow rate m. 

[*] This flow geometry is significantly different from previous cases that assume spherical 
symmetry. For magnetically controlled outflows, only a fraction of the field lines (and hence 
streamlines) allow outflow. In addition, passage through the sonic point depends sensitively 
on the divergence operator (at the transition) which depends on the configuration of the 
streamlines. As the magnetic field strength of the background (from the star) increases 
relative to that of the planet, more of the surface has open field lines, but not all of the 
streamlines allow for smooth sonic transitions. Another key difference is that spherical flows 
are often taken to be "self-limiting", where the outflow rate, in part, determines the optical 
depth to the incoming photons. In this geometry, most of the heating photons impinge upon 
the system from the equatorial directions, whereas most of the outflow is directed along the 
poles of the planet. This flow geometry also changes the effect of tidal forces exerted on the 
planet by the star: Because mass loss occurs along the polar directions, the tidal forces act 
to inhibit, rather than enhance, the mass outflow rates (see also TAL; compare with MCM). 



6.2. Discussion 

Since this problem contains a number of physical parameters, it is useful to summarize 
them here and discuss which ones are the most important. The structure of the magnetic 
field requires the specification of three quantities: the surface field strength Bp, the effective 
planetary radius Rp, and the ratio (5 of the background field strength to the surface field 
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strength. For the cases under consideration, however, the magnetic field is assumed to be 
strong enough to guide the flow, so that only the field geometry plays a role. In this simpli- 
fied treatment, the magnetic field geometry is characterized by the single parameter /3 (see 
equation [11]) . For purposes of launching the wind from the planet, the main contribution 
of the stellar field is to provide a nearly constant background field, taken here to lie in the 
vertical direction. Since the planetary field is much stronger in the region within the sonic 
surface, this assumption is valid for determining the launch of the wind. However, the stellar 
field direction will not necessarily line up with the pole of the planet, so that more general 
geometries should be considered in the future. In addition, the stellar field configuration 
will affect the manner in which the outflow propagates after its launch — beyond the sonic 
surface — and this problem should also be addressed in future work. 

With the flow geometry set by the magnetic field structure, the outflow problem deter- 
mines the fluid variables as a function of the (single) coordinate p which follows the field lines 
(Section HJ). These quantities include the density p(p), flow speed v(p), temperature T(jp), 
and ionization fraction X + (p) = n + /no, and depend on the planet mass Mp and radius Rp, 
and the stellar heating flux Fjjy. Under the assumption of isothermal flow, the problem is 
reduced further to two dimensionless fluid variables a = pj p\ and u = v/a s , along with the 
dimensionless depth of the planetary gravitational well b = GMp/(a 2 s Rp). The dimensionless 
problem thus has only two parameters (b,/3), and they determine the dimensionless outflow 
rate (see equations |57] - |65] and Figures |5] and [6]) . 

The determination of physical quantities requires specification of the density scale pi, 
the sound speed a s , and the planet radius Rp. The UV flux Fuv from the star determines, 
in part, the sound speed a s and the density p\ at the base of the flow. Over the expected 
parameter space, the sound speed a s and density scale n\ vary by only factors of 3 - 10. 
However, the outflow rate decreases exponentially with the depth b of the potential well, and 
shows complicated dependence on the magnetic field strength parameter (3 (see Figures \5\ El 
and [7]). Note that this exponential sensitivity to the depth of the potential well is a common 
feature in outflow problems where the sound speed is less than the escape speed (compare 
with the case of outflows from circumstellar disks driven by external FUV radiation; Adams 
et al. 2004). Since the dimensionless outflow rate displays exponential dependence on b, 
where b = GM P /Ra 2 , the exact value of the temperature (or equivalently the sound speed) 
can be important. In addition, the radius R in this expression corresponds to the radius R\ 
where the outflow is launched, i.e., the r = 1 surface. Although the the radius R\ differs 
from the planet radius Rp by only ~ 10 percent, this difference can be significant for the 
regime of large b where the outflow rates are exponentially suppressed. 
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6.3. Future Work 

This paper represents only the first step toward understanding outflows from Hot 
Jupiters in the regime where magnetic fields dominate the flow geometry. This work should 
be carried forward in a number of directions: 

One approximation used here is the assumption of isothermal flow. For the next stage 
of development, an analogous calculation can be carried out using a more general, polytropic 
equation of state. However, a full treatment of the heating and cooling should be undertaken. 
This calculation requires one to solve the energy equation ( 132|) and the ionization equation 
f !36p . where the heating from the central star is determined through a three-dimensional 
radiative transfer calculation. The chemistry of the outflow should also be included, both to 
get a better description of the heating and cooling processes and to determine observational 
signatures (e.g., Garcia Munoz 2007). Note that one of the intrinsic complications that arises 
in this problem is that the flow configuration (determined by magnetic field lines) does not 
have the same geometry as the heating and cooling processes. 

This paper considers a dipole magnetic field on the planet and a constant background 
field from the star; in addition, this work specializes to the case where the stellar contribution 
to the magnetic field near the planet is purely vertical. A wide range of magnetic field 
configurations are possible in star/planet systems and these possibilities should be explored 
further. One particular issue arises with sufficiently strong stellar fields, which provide a 
background field for purposes of launching the wind from the planet. If the stellar field is 
sufficiently strong and straight, the flow cannot pass throughly through a sonic transition. 
In this case, the sonic point is effectively removed to spatial infinity and the outflow rates 
are suppressed (e.g., see equation [55]). For the geometry considered here this suppression 
arises when the field strength parameter 0> 0.01 (which requires the stellar surface field to 
be ~ 10 times that of the planet for a = 0.05 AU, or equal to that of the planet for a = 0.023 
AU). A more detailed study of this regime should be undertaken, including more complicated 
configurations for the background (stellar) field. Another related effect is that the interaction 
between the planetary magnetic field and the stellar magnetosphere, including the flow of 
material considered here, can lead to an enhancement in stellar activity (e.g., Cuntz et al. 
2000, Cohen et al. 2009) and orbital evolution (Chang et al. 2010). This paper also focuses 
on the case where the stellar field is vertical, which would arise, e.g., from a dipole field from 
the star that rotates with the same angular velocity as the planetary orbit. Future work 
should relax this assumption. In particular, the field lines could trail the planet, wrap up, 
and lead to further complications. 

For sufficiently strong magnetic fields and low outflow rates, the ram pressure of the 
flow is not strong enough to affect the underlying field structure. For weaker fields and/or 
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higher outflow rates, however, the back-reaction of the flow on the magnetic field should be 
taken into account. An understanding of this physics can be attained through the solution 
to the Grad-Shafranov equation (e.g., Shafranov 1966), which self-consistently determines 
the distribution of streamlines (see also Cai et al. 2008). Even for the strong field limit, 
considered here, this back-reaction plays a role near the X-point. On a related note, the 
closed field lines can also be loaded with mass and the planet will thus develop a quasi-static 
exosphere (see TAL); this exosphere can enhance the observational signature, by increasing 
the transit depth at UV wavelengths, and should thus be considered further. 

Finally, we note that some of the planets observed in transit have orbits with nonzero 
eccentricity. The canonical example is the planet HD17156b (Barbieri et al. 2007), which 
has an orbital eccentricity of e 0.67; in this system, the stellar flux varies by a factor 
of ~ 25 over the course of the planetary orbit. The planetary outflow in this system, and 
others with similar architecture, will thus be time dependent. The resulting time dependent 
outflow rates will be sensitive functions of the heating and cooling mechanisms, as well as the 
(complicated and time dependent) magnetic field configurations. The study of such systems 
will provide sensitive tests of the outflow mechanism. 
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A. Equivalent Spherical Problem 

One complication inherent in the study of magnetically controlled outflows is that the 
problem is three dimensional. Even the reduced problem considered here, with a constant 
stellar field, requires two dimensions. In order to gain further insight into the problem, this 
Appendix develops an equivalent spherically symmetric version of the wind problem. 

The geometry of the flow pattern is specified by the divergence operator. The quantity 
that appears in the conditions that match solutions at the sonic point can be written as 



Go P ^(§)~\ (Al) 

where the geometrical factor Q is given by equation (|46p . and dl;/dp is given by equation 
fj47|) . In the limit of large £ ^> 1, one can show that G op ~ 12/(/3£ 4 ). In the opposite 
limit where £ ~ 1, near the planetary surface, G op ~ 3/£. This latter form is consistent 
with the result for a dipole field configuration near the poles of the system (recall that the 
outflow is concentrated near the poles). These two limiting forms can be connected through 
intermediate values by adopting the form 

3 

Gop = W+WM ' (A2) 

The integrated form of the dimensionless continuity equation then becomes 

au TT^T4 = x > (A3) 

where A is a (single) constant for the equivalent spherical problem. The differential form 
of the continuity equation is obtained by taking the derivative of equation fl A3[) . Near the 
planet, we thus obtain d{au^)/d^ = 0, the form appropriate for a dipole divergence. Far 
from the planet, the continuity equation reduces to the form d{au)/d^ = 0, which is the 
form applicable to flow along a single (cartesian) direction. The force equation remains the 
same (see equation 



With this form for the divergence, the condition for flow passing smoothly through the 
sonic point takes the form 

b =TTwrv < A4 > 

which thus specifies the radius £ s of the sonic transition (compare with equation [52]). In 
the limit j3 — > 0, this matching condition becomes b = 3£, as expected for a dipole. This 
expression flA4|) has no real solutions for sufficiently large b and/or large 0. The condition 
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required for solutions to exist, and hence for the flow to pass through the sonic point, can 
be written in the form 

/36 3 < 16 . (A5) 

For cases where sonic transitions are possible, the required value of the constant A is thus 
given implicitly by the relation 

l -\ 2 (l + /3/4) 2 - In A = ~ + ln(/3/4 + C 3 ) + 6(1 - V&) , (A6) 

where £ s is given by the solution to cubic equation (1A4|) . For consistency, the area subtended 
by the outflow must be proportional to the function £ 3 /(l + /3£ 3 /4) that appears in the 
continuity equation ( 1A3I) . This quantity starts near unity at the planet surface, grows like 
£ 3 near the planet, and then approaches a constant value in the limit £ — » oo. This behavior 
is thus analogous to that of the more physical problem considered in the text. 
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